Influence of humidity on granular packings with moving waUs 
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A significant dependence on the relative humidity H for the apparent mass (Mapp) measured at the 
bottom of a granular packing inside a vertical tube in relative motion is demonstrated experimentally. 
While the predictions of Janssen's model are verified for all values of H investigated (25% < H < 
80%), Mapp increases with time towards a limiting value at high relative humidities {H > 60%) 
but remains constant at lower ones (H = 25%). The corresponding Janssen length A is nearly 
independent of the tube velocity for H > 60% but decreases markedly for H = 25%. Other 
differences are observed on the motion of individual beads in the packing: for H = 25%, they are 
almost motionless while the mean particle fraction of the packing remains constant; for H > 60% 
the bead motion is much more significant and the mean particle fraction decreases. The dependence 
of these results on the bead diameter and their interpretation in terms of the influence of capillary 
forces are discussed. 

PACS numbers: 45.70.-n, 68.08.Bc,81.05.Rm 

Dense granular flows in vertical channels are encountered in many industrial processes and often display inter- 
mittency or blockage effects representing important practical problems 0, S • These effects depend on the force 
distribution in the moving grain packing and on its interaction with the walls: these depend in turn on the relative 
humidity H of the atmosphere. In static packings, capillary forces strongly influence, for instance, the stability of 
sand piles as a function of H S 0? S iS? .9] • the present work, the influence of humidity on force transmission in a 
granular packing inside a vertical tube in relative motion is analyzed from variations of its apparent weight measured 
at the bottom of the packing. 

Stress transmission in static or quasistatic granular packings has been frequently investigated theoretically and 
experimentally ^3 ] following the pioneering work and Janssen . This latter study 

and subsequent ones predict that the vertical stress at the bottom of a grain packing inside a vertical cylindrical tube 
reaches exponentially a limiting value as the height of the packing increases. This is due to the redirection of vertical 
stresses towards the side walls where friction between grains and the walls occurs: this shields the lower grain layers 
from the weight of the upper ones so that the apparent mass Mapp measured at the bottom is lower than the total 
mass M of the grains. 

In a recent paper ^3 we demonstrated that, for a constant relative humidity H = 50%, Janssen's predictions 
remain valid for a grain packing inside a vertical tube moving upwards relative to the grains at velocities between 
10~^ and 35 mm . The present study uses these results as a tool to analyse the dependence of stress transmission 
on the relative humidity H in the range 25% < H < 80%: H will be shown to influence strongly the variations and 
value of the apparent mass of the grains during the tube motion. Analyzing the motion of the grains during the tube 
displacement and the influence of the bead diameter confirms this influence. 

The experimental setup includes a vertical glass tube (length = 400 mm, internal diameter D = 30 mm) containing a 
mass M of monodisperse glass beads (diameter density p = 2.53 x 10^ kg m~^). The tube can be moved at velocities 
V between 10~^ to 130 mm over distances of 110 mm with accelerations up to 0.7 ms"^. The acceleration and 
deceleration distances are of order 12 mm at the fastest velocity and of less than 10~^ mm for V < 10 mm s~^. A 
static cylindrical piston (diameter 29 mm) is inserted at the bottom of the tube without any mechanical contact 
(Fig. it is screwed onto a sensor measuring the force applied by the packing on the piston. The apparent mass 
of the grains is determined up to 256 times per second with an uncertainty of ±0.1 g. The relative humidity H is 
regulated to within ±2%. 

Figure 121 displays variations of the apparent mass Mapp as a function of time for three different relative humidities 
H = 25%, 60% and 80% and a same diameter d = 2 mm of the beads. The total mass of grains is M = 300 g in all 
experiments (the packing height is :^ 270 mm). As expected from Janssen's model, the apparent mass after filling 
the tube [domain (1) iri Fig.[2j is lower than M. Its value Minit varies strongly from one experiment to another (typ. 
±50% of the mean value) due to the variability of the force chains redirecting the weight towards the side walls. 
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FIG. 1: Schematic view of the experimental set-up. 



As the tube starts to move upwards, Mapp decreases suddenly. In contrast with the value of Minit^he variation of 
Mapp during the motion [domain (2) in Fig.I2j is very reproducible from one experiment to another 17] at all relative 
humidities H (typical deviation ±5%). These variations of Mapp with time depend however on H (Fig.[2|). 

For H = 60%, Mapp increases towards a limiting value M^^^ while the tube moves. M^^^ is reached over a distance 
independent of the tube velocity V. This result coincides with that reported previously for H = 50% JT^. The curves 
are similar for H = 80% (Fig. ^ and, at a given V, variations for H = 60% and 80% overlay well (note that the 
difference between the values of Minit for H = 60% and 80% are not meaningful and correspond to the fluctuations 
pointed above). For H = 25%, Mapp still drops abruptly after the tube starts to move but remains practically constant 
thereafter. This constant value is also larger than that of Mf^p for = 60% and 80%. 

After the tube has stopped [domain (3) in Fig. [21, Mapp increases again for all relative humidities H: the final 
apparent mass Mfinai remains however lower than the initial one Minit and the dispersion of its values is significantly 
larger than in domain (2). 

These results indicate that lowering H influences the dynamics of the rearrangement of the force chains in the 
packing when the walls are set in motion. In order to study quantitatively these effects, we checked that Janssen's 
predictions were verified for all H values. This was realized by studying variations of the limiting value M^^p as a 
function of the total mass M. The model assumes that the redirection of the vertical stress azz into an horizontal 
component cr^r is characterized by a coefficient K and that the ratio between the tangential and normal stresses at the 
walls is the Coulomb friction coefficient ja. Assuming a zero vertical stress dzz at the surface of the packing (z = 0), 
its value crzz{zo) at the bottom {z = zq) verifies: 

azz{zo)=pcgX(l-e-'^/^y (1) 

in which the characteristic length A is given by: 

\ = D/4.iiK. (2) 




FIG. 2: Variation of the apparent mass Mapp of a grain packing {d — 2mm) of total mass M — 300^ as a function of time: 
(i) right after filling the tube, {2) tube moving upwards {V — 4mm s~^), (3) tube stopped. The curves correspond to three 
different relative humidities: (— ) H = 25%, {- -) H = 60%, {- - -) H = 80%. 
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FIG. 3: Reduced apparent mass Mfpp/Moo vs. reduced total mass M/Moo (bead diameter d = 2 mm). The solid line 
corresponds to Eq. © — (a) Tube velocity: V = 4mms"\ relative humidity: (•) H = 25%, (■) H = 60%, (a) H = 80%. 
— (b) Relative humidity: H = 60%, tube velocity: (♦) V = 0.02mms"\ (■) V = 0.2mms"\ (•) V = 2mms"\ (a) 1/ = 
20mms~\ (t) V = 34mms~^. 



The eflPective weight Mg^^^g measured by the sensor is equal to azz{zo)7rD'^ /A and the total mass M of the grains is: 
PczqttD'^ /4 (c is the solid volume fraction of the packing). By combining Eqs. (0 and ((2j, one obtains the following 
expression of the apparent mass Mappi 

where Moo is the predicted asymptotic value of Mapp when M increases: 

Moo = pcXirDyA. (4) 

In order to verify the validity of relation the ratio M^^^/M^^ is plotted in Figs. |3Ji-b as a function of the 
reduced total mass M/M^^ {M^^ is taken equal to the limiting value of M^^^ for large values of M). Both for 
different relative humidities H (Fig.[2K) and different tube velocities V (Fig.lHt)), the variations of H coincide precisely 
with theoretical predictions from Eq. 

Following these results, Janssen's hypothesis can be assumed to be valid in all our experiments. For given velocity, 
relative humidity and bead diameter values, a single experiment using a large total mass M is required to determine 
the length A: A can indeed be computed from Eq. Q using the asymptotic apparent mass M^^ measured in this 
way. The relative error on the value of A is of the order of ±4%. 

The influence of the relative humidity is then characterized by comparing the variations of A with the tube velocity 
for three different values of H (Fig.EJi). For H = 60% and 80%, A is nearly constant and close to the tube diameter 
D with A = (32 ± 1.5) mm = I.ID over a range of velocities of three decades. On the contrary, for H = 25%, A 
decreases roughly linearly with V from 65 to 35 mm over the same range of velocities (in this curve, data points for 
V 50 mm could not be obtained due to a mechanical resonance of the set-up). 

Another important parameter is the diameter d of the beads: variations of A as a function of V are displayed in 
Fig.2t> for two different relative humidities {H = 25% and H = 80%) and four different diameters d from 1.5 to Amm. 
At a given relative humidity similar variations of A with the tube velocity V are observed for all bead diameters. 
For H = 80%, A is constant or varies slowly with V, while it decreases by a factor of 2 for = 25%. However, at 
H = 25%, the variation is slower for the smallest beads {d = 1.5 mm) than for the others. 

Variations of the normalized apparent mass Mapp/M^^p with the normalized displacement z/D provide additional 
information. Overall, the same qualitative differences between H = 80% and 25% are observed for all bead diameters 
as in Fig. 121 for d = 2 mm. For H = 80%, after a sharp initial drop, Mapp increases slowly towards a limiting value 
at nearly all velocities, while it remains roughly constant for H = 25%. For the intermediate value H = 60%, the 
variation may depend however on the bead diameter. In Fig.jSl for d = 1.5 mm and d = 2 mm, Mapp/M^^p increases 
with z/D towards 1. For the largest beads {d = Amm), on the contrary, Mapp/M^^p remains roughly constant while, 
for d = 3 mm, the variation is intermediate. These observations may indicate that, for H = 60%, small beads are 
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FIG. 4: Variations of Janssen's length A with the tube velocity V — (a) for relative humidities: (•) H = 25%, (■) H = 60%, 
(a) H = 80% and bead diameter d = 2 mm. — (b) for relative humidities H = 25% (black symbols) and H = 80% (open 
symbols) and for several beads diameters: (V,T) d — 1.5mm, (o,«) d — 2mm, (□,■) d — 3mm, (A, A) d — Amm. 
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FIG. 5: Variation of Mapp/Mf^p with z/ D for a tube velocity V — A mm s ^, a relative humidity H — 60% and bead diameters: 
( ) d = 1.5 mm, ( ) d = 2 mm, ( ) d = 3 mm, (• • •) d = 4mm (M = 300^). 



initially dragged upwards by the motion of the tube and that rearrangements of the packing occur later. For larger, 
heavier beads, this effect is smaller. 

Additional clues are provided by video recordings and spatio-temporal diagrams of the bead motion close to the wall. 
At low relative humidities, they are almost motionless while, at large ones, they move permanently. The amplitude 
of the motion is larger for smaller beads and at high velocities and, in the transition regimes, in the upper part of the 
tube. 

Spatio-temporal diagrams of light intensity variations on the outside of the packing (Fig. ^ are realized with a 
linear CCD camera (sampling rate: 2000 lines/s). The x-axis corresponds to time and the 2:-axis to distance along a 
vertical line in the packing. Figure El corresponds to an experiment during which significant bead motions occurred 
(note that the acceleration of the tube is instantaneous at the scale of the diagram) . Horizontal stripes at the left and 
right ends of the diagram are obtained before and after the motion. At the beginning of the displacement, the upper 
surface of the packing rises sharply (strongly inclined stripes in the diagram) and the mean particle fraction decreases 
by Ac :^ 1.5%. After the height of the packing has become constant, the grains still drift upwards, but slower and 
a compensating downwards flow takes place in the central part of the packing. Finally, after the tube has stopped, 
the packing surface does not move down to its original location. No clear-cut dependence of these particle fraction 
variations on the tube velocity was observed. 

The upper curve in Fig. El displays the variation of Mapp with time during that same experiment: the variations 
of Mapp lasts a little longer than those of the mean particle fraction. This implies that rearrangements of the force 
distribution still occur after the decompaction has stopped. At low relative humidities {H = 25%), on the contrary. 
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FIG. 6: Spatio-temporal diagram of the top of the grain packing [z — distance from bottom of the packing, slope of dashed line 
corresponds to a 25mm velocity). Upper curve: simultaneous recording of Mapp as a function of time {V = 40 mm s~^, 
d 2mm, H 60%). 

no upwards drifts of the beads or of the packing surface are observed. 

The results presented above demonstrate that dynamical variations of the apparent mass Mapp are very sensitive 
to the relative humidity H although Janssen's model still remains valid. This dependence is likely to correspond to 
variations of the capillary forces between grains and between the grains and the walls. Many recent studies have dealt 
with the dependence of these forces on the content and distribution of liquids in granular materials j?], @, S |2l| . Even 
for small amounts of liquid at the contact point between beads, large capillary forces appear on a single asperity or 
on several asperities without requiring that pendular liquid bridges build up over the full contact region. 

At large relative humidities (typ. H = 60%, 80%), interactions between grains and with the walls are strong enough 
so that grains get dragged upwards by the tube walls (reducing strongly Mapp in the process) and keep moving even 
after a constant value of Mapp has been reached. A significant decompaction of the packing occurs (by 2.5% in our 
experiments) and does not disappear after the tube has stopped. Very similar results are obtained for H = 60% and 
80%: this implies that interaction forces do not depend very much on H once liquid bridges have built up on the 
rugosities. 

For low relative humidities {H = 25%), the capillary interactions should be weaker and take place on a few asperities: 
this would explain why beads remain almost motionless during the experiment although the interaction with the walls 
is strong enough to modify the force distribution inside the packing at the onset of the motion. Also, for single beads, 
the ratio of the capillary forces and the weight of the grains is likely to be larger for smaller beads, explaining the 
distinct variation for d = 1.5 mm in Fig. EId. The smaller amplitude of the displacements of the beads allows in 
addition to reach faster a dynamical equilibrium configuration: this is shown by the time dependence of Mapp which 
becomes immediately constant in contrast with higher values of H. Also, for H = 25%, the perturbation induced by 
the motion increases with the final velocity V of the tube: this is shown by the linear decrease of M^^^ with V and 
contrasts with its constant value for H = 60% and 80%. 

Note that no electrostatic effects were observed, even for H = 25%. 

To conclude, these experiments have allowed for a sensitive study of the dependence of force screening effects on 
the relative humidity. Increasing the latter enhances interactions between grains and with the walls, resulting in a 
significant decompaction and in a much larger amplitude of the grain motions. Janssen's model however still describes 
correctly in all cases variations of the effective mass with the height of the packing: the main difference is a decrease 
of the screening length with an increasing tube velocity for the lowest relative humidity while it is constant at larger 
ones. In future works, this information will need to be complemented by the influence of two additional important 
parameters: the roughness of the tube walls and of the beads and the mean particle fraction of the packing. 

We thank Ch. Frenois for the realization and development of the motion control system, G. Chauvin and R. Pidoux 
for the realization of the experimental setup, E. Clement and B. Perrin for helpful discussions and Ph. Gondret for a 
careful reading of the manuscript. 
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